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Hair follicles (HFs) undergo cyclical periods of growth, which are fueled by stem cells (SCs) at the base of the 
resting follicle. HF-SC formation occurs during HF development and requires transcription factor SOX9. Whether 
and how SOX9 functions in HP-SC maintenance remain unknown. By conditionally targeting Sox9 in adult 
HF-SCs, we show that SOX9 is essential for maintaining them. SOX9-deficient HF-SCs still transition from 
quiescence to proliferation and launch the subsequent hair cycle. However, once activated, bulge HF-SCs begin to 
differentiate into epidermal cells, which naturally lack SOX9. In addition, as HF-SC numbers dwindle, outer root 
sheath production is not sustained, and HF downgrowth arrests prematurely. Probing the mechanism, we used 
RNA sequencing (RNA-seq) to identify SOX9-dependent transcriptional changes and chromatin immunoprecip- 
itation (ChIP) and deep sequencing (ChlP-seq) to identify SOX9-bound genes in HF-SCs. Intriguingly, a large 
cohort of SOX9-sensitive targets encode extracellular factors, most notably enhancers of Activin/pSMAD2 
signaling. Moreover, compromising Activin signaling recapitulates SOX9-dependent defects, and Activin partially 
rescues them. Overall, our findings reveal roles for SOX9 in regulating adult HF-SC maintenance and suppressing 
epidermal differentiation in the niche. In addition, our studies expose a role for SCs in coordinating their own 
behavior in part through non-cell-autonomous signaling within the niche. 
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The renewal of interfollicular epidermis (IFE) is a constant 
process fueled by proliferation of cells in the basal 
epithelial layer. Epidermal differentiation can be easily 
followed by the unique spatiotemporal expression pattern 
of keratins and other structural proteins. Like all skin 
epithelial cells that maintain contact with basement 
membrane, basal undifferentiated keratinocytes express 
keratin 5 (K5), K14, and hemidesmosomal integrins a6 
and p4 (Fuchs 2007). When basal cells start to differenti- 
ate, they lose contact with basement membrane and 
move into suprabasal spinous layers, a process associated 
with down-regulation of basal markers and induction of 
Kl and KIO (Fuchs and Green 1980). As cells continue 
upward, granular layer cells begin to express filaggrin and 
loricrin, components of terminally differentiated, outer 
stratum corneum layers that provide the skin barrier 
(Fuchs 2007). 

In contrast, hair foUicles (HFs) imdergo synchronous, 
cycHcal rotmds of tissue growth (anagen), degeneration 
(catagen), and rest (telogen). During telogen, resident HF 
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Stem cells (SCs) can remain quiescent for weeks, enabling 
them to retain nucleotides and histones incorporated 
when in an active state. Such pulse-chase experiments 
traced the residence of HF-SCs to a niche known as the 
"bulge" (Cotsarelis et al. 1990; Tumbar et al. 2004). HF- 
SCs are typified by a group of transcription factors, in- 
cluding SOX9, NFATcl, TCF3/4, LHX2, and TBXl, which 
are not expressed in IFE or in the short-Hved, transit- 
ampUfying HF-SC progeny, called the matrix (Blanpain 
et al. 2004; Morris et al. 2004; Tumbar et al. 2004). 

Upon activation, HF-SCs at the base of the bulge — a 
region now referred to as the hair germ (HG) — divide 
rapidly and generate a pool of matrix TA cells. Several 
days after HG activation, HF-SCs in the bulge proliferate 
briefly to self-renew and form a trail of cells, which forms 
the upper outer root sheath (ORS) of HFs (Greco et al. 
2009; Rompolas et al. 2013). As the matrix and ORS grow 
dovmward, the dermal papilla (DP) is pushed fvirther away 
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from the bulge and upper ORS, prompting HF-SCs to 
return to quiescence. Matrix cells envelop and maintain 
contact with DP, enabling them to proliferate several 
times before they differentiate upward to form the hair 
and its channel or inner root sheath (IRS). The lower ORS 
remains proliferative and fuels replenishment of the 
matrix. As the destructive phase ensues and the matrix 
and lower ORS become apoptotic, spared ORS cells form 
a new SC niche: Upper ORS cells that divided the fewest 
times become the new bulge, those with slightly more 
activity become the new HG, and those that proliferated 
the most differentiate to form an inner bulge layer of 
terminally differentiated cells that secrete inhibitory 
BMP6 and FGF18 signals (Hsu et al. 2011). This keeps 
the HF in telogen until sufficient activating signals accu- 
mulate to again tip the balance and start a new hair cycle. 

Bulge SCs, ORS cells, and HG express K5/K14 and 
a6p4, but only bulge SCs express CD34 (Trempus et al. 
2003). Molecularly, HG is more similar to bulge than 
matrix and is thought to consist of "primed," short-term 
HF-SCs (Greco et al. 2009). Matrix cells exhibit reduced 
expression of HF-SC markers. Upon terminal differenti- 
ation, each lineage induces expression of distinct kera- 
tins. Most notable are the group of hair keratins expressed 
by the cortex and cuticle of the hair shaft and K6 
expressed by the hair shaft medulla and the companion 
layer that separates ORS from IRS. 

SOX9 is intriguing in that it regulates the formation of 
many cell types, tissues, and organs, including not only 
HFs but also chondrocytes, testis, heart, lung, pancreas, 
bile duct, retina, sweat gland, and central nervous system 
(Vidal et al. 2005; Nowak et al. 2008; Lu et al. 2012; for 
review, see Lefebvre et al. 2007; Sarkar and Hochedlinger 
2013). When Sox9 is conditionally ablated in embryonic 
skin, epidermal differentiation in the IFF and upper HF 
proceed unchecked, consistent with their natural lack of 
this factor, but HF-SCs do not form, and hair growth is 
arrested (Vidal et al. 2005; Nowak et al. 2008). Since no 
HF-SCs are generated in Sox9 mutant mice, it is not 



known whether and how this transcription factor regu- 
lates the behavior of HF-SCs once they are established. 

In this study, we explore SOXP's function in adult HF- 
SCs and provide compelling evidence that SOX9 is 
essential for their maintenance but does not play a critical 
role in HF lineage differentiation or affect HF-SC activa- 
tion at telogen —> anagen. Rather, upon activation, SOX9- 
deficient HF-SCs lose HF stemness and differentiate 
along the epidermal lineage. This leads to premature 
arrest of HF downgrowth, which is dependent on bulge 
and upper ORS (normally SOX9-positive) but not TA 
matrix and hair shaft/IRS (normally SOX9-negative) lin- 
eages. Using genome-wide RNA sequencing (RNA-seq) 
profiling of purified wild-type and SOX9-deficient HF-SCs 
coupled with immunofluorescence and biochemical anal- 
yses, we further show that SOX9-deficient bulge cells lose 
HF-SC characteristics. Using SOX9 chromatin immuno- 
precipitation (ChIP) in vivo and deep sequencing (ChlP- 
seq) on adult HF-SCs purified directly from their native 
niche, we show that SOX9's targets include genes encod- 
ing key HF-SC transcription factors as well as extracel- 
lular factors. Activin signaling genes are particularly 
affected, and gene ablation and our rescue experiments 
show that SOX9 functions at least partially through 
regulating this pathway. Overall, our results reveal new 
functions and mechanisms for SOX9 and suggest that 
through Activin signaling, HF-SCs are able to suppress 
epidermal differentiation, sustain ORS production, and 
return to quiescence in early anagen. 

Results 

SOX9-expiessing bulge HF-SCs contiibute 
to all HF lineages 

In telogen, SOX9 is detected in CD34-' HF-SCs, in HG 
cells, and in the terminally differentiated K6"^ inner layer 
of the bulge (Fig. lA; Supplemental Fig. lA). After anagen 
induction, SOX9 expression in the downgrowing HF 




Figure 1. Loss of SOX9 in adult HF-SCs causes 
deficient hair coat recovery. [A] Immunofluores- 
cence of SOX9 and CD34 in the telogen HFs. (B) 
Strategy to laiock out Sox9 in adult HF-SCs. (C) YFP 
and SOX9 immunofluorescence in telogen cKO HFs 
2 wk after targeting. SOX9 is gone in HF-SCs, but 
K6* inner bulge cells and YFP" cells in the CD34 
layer (arrows) still express SOX9. (D) Hair coat 
recovery of Sox9 wild-type (WT), Het, and cKO mice 
after RU486 treatment at the beginning of second 
telogen (postnatal days 60-74 [P60-P74]). Untreated 
Sox9"'f' X K15-CiePGR X R26YFP mice were also 
used as a control. Bars, 30 (im. The white dashed 
line denotes the epidermal-dermal border, and the 
blue channel is DAPI staining. (Ana) Anagen, (Bu) 
bulge, (Telo) telogen. 



GENES & DEVELOPMENT 329 



Kadaja et al. 



becomes restricted to early bulge progeny of upper ORS, 
known to maintain HF-SC stemness (Supplemental Figs. 
IB, 2B). Indeed, when subjected to lineage tracing with 
Sox9-CreER and Rosa26-fl-stop-flYPP (R26YFP) mice, 
SOX9-expressing YFP"^ cells gave rise to all adult hair 
lineages (Supplemental Fig. IC). Prior to tamoxifen, no 
YFP* cells were detected, documenting the specificity of 
Cre activation. 

Equally important, YFP* HF-SCs persisted long-term 
and were still fueling new hair regeneration 8 mo later. 
These results extended prior SOX9 studies on developing 
skin (Vidal et al. 2005; Nowak et al. 2008) and confirmed 
that SOX9* cells within the adult bulge display the 
characteristics of HF-SCs. 

Sox9 ablation in adult HF-SCs compromises 
hail coat legeneiation 

To address whether SOX9 is essential for SC maintenance 
and/or hair cycling, we bypassed SOX9's essential role in 
HF-SC establishment and conditionally targeted Sox9 
ablation in estabUshed adult HF-SCs. For this purpose, 
we used the KeratinlS {K15] promoter to selectively drive 
expression of CiePGR in HF-SCs within the bulge and 
HQ (Morris et al. 2004). By breeding to 80x9^^'^^ X R26YPP 
mice, Sox9 ablation could be monitored by YFP (Fig. IB). 
We administered RU486 in the extended second telogen 
to activate CiePGR and specifically ablate Sox9 in bulge 
HF-SCs (CD34-') and HQ cells (CD34") (Fig. IC; Supple- 
mental Fig. 2A,B). YFP" bulge cells included a few CD34+, 
K15* HF-SCs but were mostly terminally differentiated 
K6-I-, K15~ inner bulge cells (Hsu et al. 2011). Based on 
quantitative RNA PCR (qPCR) of mRNA from HF-SCs 
purified by fluorescence-activated cell sorting (FACS), 
Sox9 levels were decreased ~80%-95% in YFP-" HF-SCs 
compared with heterozygous (Het) controls (Supplemen- 
tal Fig. 2C). SOX9 expression in YFP" HF-SCs remained 
high, as expected. 

Loss of SOX9 in HF-SCs did not result in an overt 
phenotype as long as HFs remained in telogen. Two 
weeks after RU486 treatment, percentages of YFP"^ 
Sox9-null HF-SCs in skin were comparable with 
those of YPP"^ Sox9 Het (control) HF-SCs. In addition, 
throughout telogen, Sox9-null HF-SCs remained qui- 
escent like control HF-SCs (Supplemental Fig. 2D,E). 
As judged by shaving the skin and monitoring skin 
darkening, HFs lacking SOX9 in their SCs were able to 
enter into anagen (Fig. ID). The transition from pink 
(telogen) to black (anagen) skin is a reflection of melanin 
production by differentiating melanocytes, which occurs 
in anagen and leads to melanin transfer to concomitantly 
differentiating hair shaft precursors (Miiller-Rover et al. 
2001). 

After HFs transitioned from telogen to anagen, Sox9 
conditional knockout (cKO) mice began to show de- 
ficient hair coat recovery with patchy and incomplete 
hair growth (Fig. ID). Untreated Sox9 fl/fl;K15-CrePGR 
R26YPP mice showed normal hair coat recovery, as 
did RU486-treated Sox9 Het and wild-type controls 
(shown in Fig. ID). 



Sox9 cKO HPs launch the haii cycle normally, 
but once activated, their HP-SCs lose stemness 

Following Sox9 ablation, mutant HF-SCs generated 
a seemingly normal early anagen-phase HF and hair shaft. 
Closer inspection, however, revealed that by Anagen IV, 
the bulge region showed signs of deformation, and HF 
dovmgrowth had markedly slowed. Moreover, by full 
anagen (Anagen VI), the hair shaft and IRS, which had 
formed normally at earlier stages, now showed signs of 
increased deformation (Fig. 2A,B; Supplemental Fig. 3A- 
C). Since the morphology of early anagen cKO HFs was 
normal and since Sox9 expression is normally silenced 
after wild-type ORS cells reach the TA matrix compart- 
ment, the later stage distortions of cKO HF/shaft mor- 
phology were likely secondary due to the physical barrier 
imposed by an increasingly aberrant bulge. Consistent 
with this notion, SOX9-deficient matrix cells were still 
proliferative (Fig. 2C) and gave rise to all HF lineages, 
including the K6"^ companion layer, hair keratin"^ (mAb 
AE13*) cells of cortex and cuticle, and trichohyalin* (mAb 
AEIS"") cells of IRS and hair shaft medulla (Fig. 2D,E). 

We did notice that cells expressing HF differentiation 
markers diminished as anagen progressed. EdU labeling 
and quantification revealed that even though percentages 
of labeled TA matrix cells were normal, their total 
numbers were also significantly reduced in the Sox9 
cKO HFs (Fig. 2F). This was not attributable to apoptosis, 
as activated CaspaseS labeling was still very low in the 
mutant matrix (Fig. 2G). Pinpointing the defect further, 
we showed that SOX9-deficient HFs entered anagen with 
kinetics comparable with their wild-type counterparts, 
indicating that HQ progenitors at the bulge base still 
seemed to be responding normally to their niche cues. 
Moreover, at Anagen III, when wild-type bulge HF-SCs 
self -renew (Hsu et al. 2011), Sox9 cKO bulge HF-SCs 
behaved normally. Surprisingly however, once activated 
to self-renew, SOX9-deficient bulge HF-SCs failed to 
return to quiescence (Fig. 2C, H). Soon thereafter, HF 
downgrowth was arrested. 

Taken together, these results were all consistent with 
the view that SOX9 is essential for maintaining the bulge 
HF-SCs and upper ORS, and, without SOX9, these nor- 
mally SOX9* cells are unable to fuel the normally SOX9" 
TA matrix and differentiating hair cells. 

Within the bulge niche, activated Sox9 cKO 
HP-SCs fail to repress epidermal differentiation 

We returned to early anagen phase, where the first 
morphological consequences of Sox9 ablation in the 
bulge had emerged, which were soon followed by distor- 
tions in normally SOX9^ hair cell populations. Closer 
inspection revealed the presence of keratinized pearls 
within the early SOX9-deficient bulge (Fig. 3A). Immtmo- 
fluorescence microscopy further revealed the presence of 
KIO and filaggrin, indicative of terminal differentiation 
(Fig. 3B,C). The aberrant IFF differentiation appeared to be 
confined to the niche, as when Sox9 cKO HFs entered 
into anagen, YFP* Sox9 cKO progeny along the ORS 
showed no signs of KIO and filaggrin immunostaining. 
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AnaVI AnalV AnaVI Anagen VI post-RU486 




Figure 2. Sox9 ablation in adult HF-SCs results in bulge deformation and an anagen-dependent failure to return to quiescence. [A] 
Semithin sections of Sox9 wild-type (WT) and cKO HFs in the first anagen phase of the hair cycle following RU486 treatment in 
previous telogen. Note that (1) cKO HFs are shorter, and (2) IRS/hair shaft perturbations are subsequent to those of the bulge. [B] Anti- 
SOX9 immunofluorescence. (C) EdU incorporation (24-h pulse) in HF-SCs in first (full) Anagen VI following RU486. (D,£) K6 in the 
companion layer (Cp), hair keratins (AE13) in the hair shaft (HS) cortex and cuticle, and trichohyalin (AE15) in the IRS and HS medulla 
are still expressed in Sox9 cKO HFs. [P] Quantifications of EdU incorporation and total matrix cell number of Anagen VI HFs post- 
RU486. N = 3 mice for each genotype. (G) Quantification of activated caspase 3 in the matrix of Anagen VI HFs post-RU486. N= 3 mice 
for each genotype. [H] Quantifications of EdU incorporation in the bulge at the indicated time points. N = 4 mice for each genotype. 
Data are mean ± SEM. (***') P < 0.0001 . Bars, 30 jjim. White dashed lines denote the epidermal-dermal border, and the blue channel is 
DAPI staining. Asterisk indicates the nonspecific staining of SOX9 Ab. (Amp) Arrector pili muscle; (Bu) bulge; (DP) dermal papilla; (IRS) 
inner root sheath; (ORS) outer root sheath; (SG) sebaceous gland; (Mx) matrix; (ns) not significant. 



Intriguingly, while the actively cycling (newest) cKO 
bulge was undergoing these gross morphological and 
molecular changes, the old cKO bulge, which is not 
reused after its initial hair cycle (Hsu et al. 2011), 
remained quiescent and appeared to be largely intact 
morphologically (Supplemental Fig. 3A). It also displayed 
the characteristic SOX9"^K6"^ iimer bulge layer, which still 
anchored to the old club hair (Fig. 3D). These findings 
further underscored the phenotypic dependency on hair 
cycle activation. 

Spontaneous epidermal differentiation first appeared in 
the telogen bulge and was marked by a few KIO"^ cells 



within the CD34'" HF-SC layer (Fig. 3E). At this early 
stage, Sox9 cKO HF-SCs still exhibited molecular signs of 
stemness (Supplemental Fig. 2A). Soon after anagen entry, 
however, CD34 and other markers of stemness began to 
wane (Fig. 3E), and by the end of the first hair cycle in 
which Sox9 had been depleted, HFs were reduced to 
epidermal cysts (Fig. 3F). 

Taken together, we placed the major SOX9 defect 
downstream from HF-SC survival and anagen onset but 
at the time (mid-anagen) when activated bulge HF-SCs 
normally return to quiescence. Moreover, our data showed 
that in SOX9's absence, activated HG progenitors still 
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commit to HF differentiation programs, but activated 
bulge HF-SCs fail to both return to quiescence and also 
suppress epidermal differentiation within the niche. 

Wnt and Notch signaling still occur in Sox9 cKO HPs 

KIO* cysts have also been described when Notch signal- 
ing is compromised in HFs (Demehri and Kopan 2009; Hu 
et al. 2010; Lin et al. 2011). However, although signs of 
active Notch signaling, including the nuclear Notch 
intracellular domain (NICD) and downstream NICD 
target gene Hesl, were detected in the lower bulge/HG 
of telogen-phase HFs, they showed no significant differ- 
ences between wild-type and cKO HFs (Supplemental Fig. 
4A,B). Thus, changes in Notch signaling did not seem to 
account for the phenotype that we observed. 
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KIO"^ cysts have been also reported in HFs expressing 
a LEFl mutant that lacks its p-catenin-binding domain 
and hence cannot be transactivated in response to Wnt 
signaling (Niemann et al. 2002). However, in the absence 
of SOX9, Wnt signaling still occurred, as judged by 
nuclear |3-catenin in the activated HG cells and LEFl in 
the emerging matrix cells (Supplemental Fig. 4C,D). 
Additionally, mRNA levels of Wnt targets seemed to be 
unaltered at the time point when phenotypic changes in 
Sox9 cKO telogen HFs were already evident (Supplemen- 
tal Fig. 4E). mRNA levels of Tcf3 and Tcf4 were reduced 
somewhat but not enough to suggest that Wnt is re- 
sponsible for the changes that we observed in the Sox9 
cKO bulge. 

Finally, inflammation can also lead to hair loss and cyst 
formation. Flow cytometry analysis, however, did not 
reveal any signs of immune cell infiltration in Sox9 cKO 
skin (Supplemental Fig. 4F). 

HF-SCs that lose SOX9 show decreased expression 
of HP-SC signature genes 

Delving deeper into how SOX9 loss leads to a HF-to-IFE 
lineage fate switch, we transcriptionally profiled late 
telogen-phase HF-SCs (YFP-'/Scar/CD347a6hi) purified 
from Sox9 cKO and Het mice treated 2 wk earlier in 
telogen with RU486 (Supplemental Fig. 5A). We chose 
this timing because typifying sternness markers were still 
expressed by SOX9 -deficient HF-SCs, but early molecular 
changes were detectable (Supplemental Fig. 2). 

mRNAs isolated from these two samples were sub- 
jected to deep sequencing and bioinformatics analyses. 
Two-hundred-seventy-five protein coding genes were 
differentially expressed by >1. 5-fold (P-value < 0.05) in 
Sox9-null versus Het HF-SCs. Interestingly, mRNAs that 
were down-regulated by >1.5X in Sox9-null HF-SCs 
overlapped significantly with not only the HF-SC signa- 
ture but also the HF lineage, while mRNAs up-regulated 
by > 1.5-fold in Sox9-null HF-SCs overlapped most signif- 
icantly with mRNAs enriched in basal IFF cells (Fig. 4A; 



Figure 3. Sox9 cKO HF-SCs differentiate into IFE cells. {A] 
Ultrastructure of the SC niche of wild-type (WT) and Sox9 cKO 
bulge of Anagen VI HFs after RU486 treatment. The new bulge 
(Bu), deformed at this stage, is marked by red pseudocolor; dark 
brown labels the old bulge, and light brown indicates the 
sebaceous gland. The boxed area, magnified at right, reveals 
keratinized pearls and keratohyalin granules (arrows) within the 
cKO bulge. (B) KIO. (C) Immunostaining for filaggrin (FIL). KIO 
reveals atypical labeling of epidermal markers in the cKO bulge. 
(D) Immunolabeling shows that the K6* inner bulge layer is not 
targeted and remains intact during the first anagen following 
Sox9 ablation in HF-SCs. Dashed lines denote the epidermal- 
dermal border. (£) Representative images of KIO immunostainings 
of telogen end-phase HFs treated with RU486 3 mo earlier at 
telogen start. Quantifications are at riglit. N = 3 mice for each 
genotype. Data are mean ± SEM. ( * * * * ) P < O.OOOI. [F] Depletion 
of SOX9 in adult HF-SCs leads eventually to the formation 
of KlO-positive cysts. Blue channel is DAPI staining. (Cntrl) 
Control. Bars, 30 ixm. 
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Figure 4. Changes in gene expression in HF- 
SCs lacking SOX9. [A] Transcriptional profile 
of SOX9-deficient HF-SCs is shifted away 
from the normal HF-SC "signature" toward 
the basal IFE "signature." The top diagram 
shows overlap of genes differentially expressed 
by al.SX in Sox9 cKO versus Het HF-SCs 
[Sox9 RNA-seq), and genes were enriched 
>2x in telogen HF-SCs versus basal IFE 
(HF-SC signature genes). The bottom dia- 
gram shows overlap of Sox9 RNA-seq and 
genes enriched >2x in IFE versus telogen 
HF-SCs (IFE signature genes). [B] Table shows 
examples of genes >1. 5-fold down-regulated 
(green) or > 1.5-fold up-regulated (red) in Sox9 
RNA-seq. (C) Gene ontology (GO) analyses 
of genes differentially expressed by >1.5X or 
>2x in Sox9 cKO versus Het HF-SCs denotes 
the "extracellular region" category as the most 
significantly enriched compared with the 7.6% 
frequency in the Mouse Genome Informatics 
(MGI) database. (D) KIO and YFP immunoflu- 
orescence of telogen Sox9 cKO bulge prior to 
hair cychng. (Blue) DAPI. At right, the per- 
centage of K10+ cells within YFP+ and YFP" 
HF-SCs of the mosaic Sox9 cKO bulge (YFP is 
in 60%-80% of HF-SCs) and wild-type (WT) 
bulge. N = 3 mice for each genotype. Data are 
mean ± SEM. (****) P < O.OOOI. Bar, 30 |j.m. 



Supplemental Fig. 6B; Supplemental Table 2). The switch 
included repression of a cohort of HF lineage transcrip- 
tional regulators, including HF-SC genes Lhx2, Sox9, 
Tcf7l2, Tbxl, and Nfatcl (Fig. 4B; Supplemental Fig. 6A; 
Supplemental Table 2). Additionally, transcriptional reg- 
ulators for HF differentiation were transcribed at even 
lower levels than normal. This shift away from HF-SC 
and toward the IFE signature was consistent with the IFE 
differentiation observed within the Sox9 cKO HF bulge. 

Gene ontology analyses revealed that many SOX9- 
dependent changes in gene expression were associated 
with secreted proteins (Fig. 4C). This was noteworthy, 
since the best-known SOX9 transcriptional targets of 
chondrocytes include extracellular matrix (ECM) genes 
such as Col2al (Lefebvre et al. 1997). To test the possi- 
bility that the SOX9-dependent phenotype may at least in 
part be rooted in cell-nonautonomous defects, we quan- 
tified KIO expression in YFP"" Sox9' and YFP" 80x9" HF- 
SC subpopulations from Sox9-cKO X R26YFP mice. As 
shown in Figure 4D, KIO was detected in both bulge 
populations irrespective of whether SOX9 was present. In 
contrast, KIO was not seen in the bulge of wild- type HFs 
(Fig. 4D, shown in the graph). Together, these findings 
suggest that within the bulge niche, Sox9-null HF-SCs 
secrete factors that induce IFE differentiation in them- 
selves and their neighboring bulge cells. 

SOX9 functions primarily as an activator 
of gene expression in adult HF-SCs 

Second telogen HF-SCs could be FACS-purified in suf- 
ficient quantities to conduct SOX9 ChlP-seq in vivo 



(Supplemental Fig. 5B). One-thousand-three-hundred- 
eighty-four genes contained SOX9-binding sites in either 
the promoter (±2 kb of transcription start sites) or 
"enhancer" (—50 kb upstream of and -i-5 kb downstream 
from the gene body) (Fig. 5A; Supplemental Table 3). 
This list included Nfatcl, Lhx2, and Tcf7l2 (TCF4), 
previously shown to be critical for maintaining HF-SCs 
(Merrill et al. 2001; Nguyen et al. 2006, 2009; Rhee et al. 
2006; Horsley et al. 2008; Folgueras et al. 2013). In this 
regard, SOX9 resembled SOX2, which, in embryonic SCs, 
binds to genes encoding pluripotency regulators (Boyer 
et al. 2005; Chen et al. 2008). 

Comparison of SOX9 ChlP-seq data with our list of 
genes > 1. 5-fold differentially expressed in Sox9 cKO 
versus wild-type HF-SCs produced a short Hst of 77 putative 
SOX9-regulated HF-SC genes (Supplemental Table 2C). 
Notably, the overlap with ChlP-seq data was significantly 
higher with genes whose expression was decreased in the 
absence of Sox9, suggesting that SOX9 functions primar- 
ily as a transcriptional activator in HF-SCs. Thus, whereas 
33% of down-regulated genes bound SOX9 (P-value < 
0.0001), only 17% of up-regulated genes did so (P-value = 
0.0004) (Fig. 5B). 

Fifty-nine percent of putative SOX9 targets that were 
down-regulated > 1. 5-fold in Sox9 cKO HF-SCs were also 
emiched in wild-type HF-SCs compared with IFE pro- 
genitors. This included Wwp2, which was previously 
identified by ChlP-qPCR as a direct downstream SOX9 
target in chondrocytes (Nakamura et al. 2011). Further 
insights into SOX9's role as a transactivator came from 
comparing our SOX9 ChlP-seq profile with activating and 
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Figure 5. SOX9 functions primarily as a transcrip- 
tional activator in adult HF-SCs. [A] Distribution of 
SOX9 peaks between promoters (—2 kb to +2 kb of 
transcription start sites) and "enhancers" (—50 kb of 
transcription start sites, gene body, -h5 kb down- 
stream from the gene) of SOX9-bound genes. (B) 
Diagram showing overlap of in vivo SOX9-bound 
targets (ChlP-seq) and mRNAs differentially 
expressed (RNA-seq) in Sox9 cKO versus Het HF- 
SCs. Overlaps represent the number of genes bound 
by SOX9 and whose expression is up-regulated or 
down-regulated by a 1.5-fold in Sox9 cKO. Overlaps 
with HF-SC signature and IFE signature genes are 
also shown. (C) ChlP-seq profiles of the S100A4 
genomic locus in wild-type (WT) telogen HF-SCs. 
Input represents the sequencing profile of genomic 
DNA used in ChIP pull-downs. Black arrows in- 
dicate locations of DNA fragments enriched in 
SOX9 Chip The SOX9 ChlP-seq profile is also 
compared with published data of histone modifica- 
tions in telogen HF-SCs (Lien et al. 201 1): H3K4me3 
Ab pulls down promoter regions, H3K79me2 marks 
actively transcribed DNA, and the H3K27me3 label 
is enriched in polycomb-repressed chromatin. Below 
is an RNA-seq profile of the S100A4 locus in Sox9 
cKO and Het HF-SCs. The locations of RNA-seq 
fragments match with exons in the S100A4 gene. All 
tracks are set to the same scale. The lengths of genes 
and direction of transcription are indicated by blue 
arrows. (D) SOX9-binding motifs (from left): the 
most enriched motif within SOX9 peaks in HF- 
SCs, the motif enriched most within SOX9 peaks 
associated with differentially expressed genes in 
Sox9 cKO versus Het, and the previously reported 
SOX9 in vitro motif (Mertin et al. 1999). Below, 
relative luciferase activities driven by SOX9 motifs 
were used as an enhancer. Assays were performed in 293FT cells {SOX9-negative) in the presence of a SOX9 expression vector or empty 
vector. All results were normalized to luciferase activity when firefly luciferase reporter with minimal CMV was used and to the 
Renilla luciferase control reporter driven by the HSV TK promoter. Experiments for each construct were performed in triplicate in 
duplicate wells. Data are mean ± SEM. (ns) Not significant; (*) P < 0.05. 
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repressive histone methylation marks on HF-SC chroma- 
tin (Supplemental Table 3; Lien et al. 2011). Figure 5C 
shows these profiles for S100A4, a representative HF-SC 
signature gene that bound SOX9. As illustrated, at least 
three distinct SOX9 peaks were identified within 
and upstream of the gene body, whereas histone H3 
trimethylation of lysine residue 4 (H3K4me3), anal- 
ogous to the H3K79me2 mark, characteristic of tran- 
scriptional elongation (Lien et al. 2011), was prevalent 
over the S100A4 gene body. Conversely, S100A4 lacked 
the repressive Ezh2/polycomb complex mark H3K27me3. 
These marks agreed well with our S100A4 RNA-seq pro- 
files, where transcription of all three exons of S100A4 was 
repressed upon loss of SOX9 in HF-SCs (Fig. 5C). 

We carried out unbiased sequence motif searches first 
within the 600 SOX9 peaks with the strongest ChlP-seq 
signals and then within peaks of genes whose expression 
changed by >1.5X upon Sox9 ablation. The resultant top 
two motifs suggested a low sequence specificity of SOX9 
binding in HF-SCs (Fig. 5D). Moreover, even the motif 
underlying the strongest peaks in HF-SC chromatin di- 



verged at several sites from the SOX9 DNA-binding motif 
that emerged from in vitro studies in which recombinant 
SOX9 was exposed to a mixture of random oligonucleo- 
tides (Mertin et al. 1999). Similar differences have been 
noted previously, first for the 44 putative targets bound by 
SOX9 in Chip hybridizations of 80 selected chondrocyte 
promoters (Oh et al. 2010) and then for the 109 putative 
targets identified in a genome-wide SOX9 ChlP-on-ChIP 
of embryonic testis chromatin (Bhandari et al. 2012). That 
said, none of these consensus sites were identical, further 
pointing to the weak specificity of SOX9 for its in vivo 
targets. Low binding specificity is emerging as a general 
feature of Sox proteins (Lefebvre et al. 1998; Kamachi 
et al. 1999). 

Chromatin binding does not necessarily imply tran- 
scription factor regulation. To address whether SOX9 has 
the capacity to regulate the SOX9-binding motifs that we 
identified from our in vivo SOX9 ChlP-seq study, we 
compared the relative activities of enhancers composed 
of 3 X multimerized motifs identified from either in vitro 
oligonucleotide analyses (Mertin et al. 1999), genome- 
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wide analyses of the strongest SOX9 peaks, or analyses of 
the SOX9 peaks from the 103 genes that exhibited 
differential expression upon Sox9 ablation in HF-SCs. 
To evaluate SOX9 dependency and enhancer activity, w^e 
transfected our expression vectors with or without an 
additional SOX9 expression vector (controlled for equal 
DNA levels). 

Each enhancer showed some SOX9-dependent activity 
(Fig. 5D). However, when compared relative to the 
Renilla luciferase internal control, the most potent en- 
hancer and the one with greatest sensitivity to SOX9 
(>15x) was that constructed from the in vivo motif that 
was present specifically in genes that also showed SOX9 
transcriptional sensitivity in HF-SCs (shovm in Fig. 5D). 
These results show that when multimerized, this SOX9- 
binding motif is sufficient on its own to act as an 
enhancer. That said, the weaker activity of the enhancer 
composed of in vitro optimized SOX9-binding sites, 
coupled with the relative lack of specificity of the in vivo 
SOX9 motifs, raises the possibility that the SOX9-de- 
pendent regulation conferred by in vivo SOX9 motifs may 
be complex and involve associated factors. Although 
beyond the scope of the present study, our findings are 
consistent with accumulating evidence suggesting that 
Sox factors do not act alone in their function (for review, 
see Kamachi and Kondoh 2013). 

SOX9 is a positive regulator of genes involved 

in Activin/TGFp signaling 

To gain further insights into the pathways preferentially 
regulated by SOX9, we focused on highly expressed HF- 
SC signature genes (FPKM >4 in Sox9 Het RNA-seq) that 
were strongly bound by SOX9 in their promoter/enhancer 
and also down-regulated in Sox9 cKO HF-SCs more than 
twofold in both RNA-seq replicates. Among this short list 
of 13 genes were Wwp2, S100A4, Sulf2, and Inhbb, 
previously reported to enhance TGFp/Activin signaling 
(Fig. 6A; Matsuura et al. 2010; Soond and Chantry 2011; 
Zheng et al. 2013). qRT-PCR on independently purified 
HF-SCs confirmed that their mRNAs were down- 
regulated in cKO HF-SCs (Fig. 6B). GAPDH and Ppib 
were internal mRNA loading controls. Bulge-specific 
reduction of S100A4 was further confirmed by immu- 
nolabeling (Fig. 6C). 

To further validate our SOX9 ChlP-seq data, we per- 
formed SOX9 ChlP-qPCR on HF-SCs (SOX9"') and IFE 
(SOX9^). These analyses confirmed that the enrichment 
obtained from SOX9 ChIP was SOX9-specific and did not 
occur when SOX9~ cells were used in this assay (Fig. 6D). 
Enrichment of DNA fragments in SOX9 ChlP-qPCR was 
specific to regions where SOX9 ChlP-seq peaks were 
detected, as shown by analyses of neighboring chromatin 
regions and random chromatin regions lacking SOX9 
peaks (shown in Fig. 6D). 

To verify that SOX9 regulates the transcription of these 
genes, we carried out luciferase reporter assays in cul- 
tured cells with or without Sox9 expression vector. 
Relative to the control minimal CMV promoter, each 
300- to 600-base-pair (bp) element encompassing SOX9- 



binding sites of either S100A4, Wwp2, Sulf2, or Inhbb 
displayed marked enhancer activity, elevating luciferase 
expression by as much as lOx in the presence of SOX9 
(Fig. 6E). The effects were SOX9-dependent, as vdthout 
SOX9, levels plummeted to that of the minimal promoter 
(Fig. 6E, red line). These data establish Wwp2, Inhbb, 
S100A4, and Sulf2 as bona fide targets of SOX9. 

SOX9 maintains the identity of HF-SCs through 
Activin signaling 

pSMAD2, the downstream effector of TGF|5/Activin 
signaling, is up-regulated at the telogen -> anagen tran- 
sition and aids in driving HF regeneration (Oshimori and 
Fuchs 2012). Since Wwp2, Inhbb, S100A4, and Sulf2 have 
been implicated as enhancers of TGFp/Activin signaling, 
their dovm-regulation in Sox9 cKO HF-SCs suggested 
that SOX9 may be involved at least in part in regulating 
this pathway. Indeed, in comparison with the wild-type 
bulge at telogen — > anagen, pSMAD2 staining was 
significantly reduced in Sox9 cKO bulge HF-SCs at the 
same stage (Fig. 7 A). 

Mice whose skin epithelium lacks TGF-(3 receptor II 
(and hence TGF(3 signaling) do not exhibit gross alter- 
ations during normal homeostasis (Guasch et al. 2007; 
Oshimori and Fuchs 2012). In contrast, when Activin 
receptor type IB [Acvrlb] is conditionally targeted for 
ablation, HFs degenerate, resulting in large cysts with 
keratinaceous debris (Qiu et al. 2011). Given these 
phenotypic parallels, we focused on the hypothesis that 
SOX9 acts on genes of the Activin pathway to sustain 
stemness and block IFE differentiation. 

Similar to Sox9 cKO HFs, cysts in Actvrlb cKO mice 
showed signs of IFE differentiation and KIO expression 
(Fig. 7B). To test whether Activin signaling is directly 
important for repressing IFE fate in the HF-SC niche, we 
asked whether Activin B, encoded by Inhbb, could rescue 
this defect. To this end, we injected RU486-treated Sox9 
cKO mice in their second telogen with recombinant 
Activin B or BSA (control) along vdth fluorescent beads 
to mark intradermal injection sites. Since prominent 
signs of HF — > IFE transition surface at telogen anagen, 
we then waxed back skins to remove hair shafts and 
induce synchronized anagen. 

As shown in Figure 7C, Activin B restored nuclear 
pSMAD2 and significantly reduced the number of KIO"^ 
iSoxP cKO HF-SCs. This decline in KIO was clearly 
specific to Activin B, as it tapered off further from the 
injection site (shown in Fig. 7C). Concomitantly, CD34 
immunofluorescence intensity was increased in Activin- 
treated, Sox9 cKO HFs. While CD34 repression is a com- 
mon sign of HF-SC proliferation, proliferation was not 
affected. 

To further study the role of Activin in SOX9-mediated 
control of the bulge niche, we treated telogen 5ox9 cKO 
HFs with Activin B and BSA, FACS-purified HF-SCs from 
injection sites, and measured their mRNA levels of 
secreted factors whose expression was altered by SOX9 
loss. Interestingly, mRNA levels of HF-SC signature 
genes Sostdcl and Sulf2 were increased, while levels of 
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Figure 6. SOX9 is a positive regulator of genes involved in Activin/TGFp signaling. [A] Top list of putative SOX9 targets that were 
highly expressed in adult HF-SCs (FPKM >5 in Het RNA-seq), strongly influenced by SOX9 loss (down-regulated by twofold or more in 
cKO versus Het], and associated with a strong SOX9 peak in ChlP-seq (peak height >13). [B] RT-qPCR verification of RNA-seq data. 
Independent FACS-purified samples were analyzed for putative SOX9 target genes involved in TGFp/Activin signaling. For RT-qPCR, 
all results were normalized to Ppib expression; GAPDH was used as an internal control. Data are mean ± SEM. N = 4 mice per 
genotype, (C) S100A4 immunofluorescence in Sox9 cKO and wild-type (WT) F-[Fs in the first full anagen after targeting. White lines 
denote the epidermal-dermal border, (Blue) DAPI, Bar, 30 (jim, (D) Independent SOX9 ChlP-qPCR on FACS-purified HF-SCs (SOX9- 
positive) and IFE (SOX9-negative) cells. For each gene, primers were designed to amplify region where the SOX9 peak was detected 
(black bars) and neighbor regions 4-8 kb away from the SOX9 peak (gray bars). Note that enrichment of DNA fragments was specific to 
HF-SCs and regions underneath the peaks, validating the specificity of SOX9 Ab and binding. Results were normalized to the intergenic 
Chr3 region; the intergenic Chr5 region was used as a control. Experiments were done in triplicate. Data are mean ± SEM, (*) P < 0.05; 
{**] P < O.OI . [E] Relative luciferase activities driven by sequences (300-600 bp) encompassing SOX9 ChlP-seq binding sites. Results are 
presented as the enhancement of luciferase activity in the presence of a SOX9 expression vector compared with an empty expression 
vector. All results are normalized to the luciferase activity of Renilla luciferase control reporter vector. Experiments for each construct 
were performed in four replicates in duplicate weUs, Data are mean ± SEM. (*'*) P < 0,001; (****) p < 0.0001, 



IFE signature genes Thbs2 and Guca2a were decreased 
after the Activin B injection (Fig. 7D). In a good correla- 
tion, Sostdcl and Sulf2 expression was decreased, while 
Thbs2 and Guca2a were increased after knocking down 
Inhbb in cukured wild-type HF-SCs (Fig. 7E). 

Once activated, Sox9 cKO HF-SCs cannot return to 
quiescence. To test whether restoring Activin signaling 
can restore this defect, we synchronized the hair cycle by 
waxing mice and then injected Activin B and BSA in- 
tradermally between 4 and 6 d after depilation; i.e., the 
time when proHferation of wild-type HF-SC starts to 
decrease. EdU labeHng and quantifications showed that 
Activin B did not affect HF-SC proHferation in either 
wild-type or cKO mice (Fig. 7F). 



Taken together, these findings were consistent with the 
non-cell-autonomous effects that we observed previously 
between SOX9 loss and IFE fate conversion and provide 
strong evidence that this link operates at least partially 
through Activin signaling. However, it also illustrates the 
multifaceted ways in which SOX9 regulates HF-SC be- 
havior. 

Discussion 

Cell-autonomous and non-cell-autonomous functions 
for SOX9 in governing the SC niche 

When Sox9 is targeted in development, HF-SCs never 
form, resulting in the arrest of SOX9-dependent HE 
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Figure 7. SOX9 maintains the identity of HF-SCs through Activin signaUng. [A] pSMAD2 immunofluorescence of HFs during the first 
telogen-anagen transition after RU486. Quantifications are at right. N = 5 mice per genotype. Data are mean ± SEM. (* *) P < 0.01. [B] 
KIO immunohistochemistry of cKO Acvrlb mice and fl/fl control (wild-type [WT]) mice. Asterisk indicates melanin. (C) Intradermally 
injected Activin B inhibits the irans-differentiation of HF-SCs to IFE caused by SOX9 loss. Mice treated with RU486 in second telogen 
were depilated to synchronize hair cycles and expedite the phenotype. Recombinant Activin B and BSA were injected intradermally 
into the same mouse at two different locations over the 3 d. Immunostaining for pSMAD2 and CD34 in Activin B and BSA-treated HFs 
are shown at left; KIO is shown at right. Quantifications of KIO* cells, EdU* cells per HF, and the total fluorescence intensity of the 
CD34 label are shown in the middle. N = 4 mice. Data are mean ± SEM. (* *) P < 0.01; (* **) P < 0.001. White lines denote the epidermal- 
dermal border. (Blue) DAPI. Bar, 30 jjim. (D) RT-qPCR quantification of mRNA levels of secreted factors in Activin B-treated bulge. HF- 
SCs were FACS-purified from Activin B and BSA injection sites independently from two Sox9 cKO mice. Data are mean ± SEM. (£) RT- 
qPCR quantification of mRNA levels of secreted factors in cultured wild-type HF-SCs, where Activin B. expression has been reduced by 
Inhbh shRNA. Data are mean ± SEM based on four independent experiments. (P) Quantifications of EdU* cells in wild-type and Sox9 
cKO bulge at Anagen IV, after treatment with BSA or Activin B. N = 2 mice. Data are mean ± SEM. 



morphogenesis but not SOX9-independent IFE (Vidal 
et al. 2005; Nowak et al. 2008). Since HF-SCs fail to 
develop without SOX9, SOX9's function as an HF-SC 
transcriptional regulator has remained unclear. This 
issue is also largely unexplored for most SOX9* cell types 
that contribute to adult tissue homeostasis (Furuyama 
et al. 2011; Guo et al. 2012; Lu et al. 2012). We now 
address this issue for the HF. Our results show clearly 
that when Sox9 is ablated in adult HF-SCs, they are 
unable to maintain stemness. Moreover, by governing 
expression of a number of extracellular factors, SOX9 can 



exert both cell-autonomous and non-cell-autonomous 
effects on SC behavior. Indeed, one of the first signs of 
Sox9 ablation was that a few intact wild-type HF-SCs 
within the niche began to aberrantly differentiate along 
the IFE lineage. 

Defects in SC behavior were manifested soon there- 
after. While HFs transitioned from telogen —> anagen and 
SOX9 -deficient bulge HF-SCs launched self-renewal on 
cue in early anagen, they failed to return to quiescence in 
Anagen IV. Despite this activity, cKO HF-SCs could not 
sustain ORS or TA matrix pools, grinding hair lineage 
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differentiation to a halt and underscoring the paucity of 
sternness within these prohferating Sox9-ablated bulge 
cells. Indeed, SOX9-deficient HF-SCs differentiated along 
the IFE lineage, fully exhausting the HF-SC pool in bulge 
and upper ORS, where SOX9 is normally expressed. 

Additionally noteworthy was the paucity of HF differ- 
entiation defects even though HF-SC-mediated fueling of 
TA matrix progenitors was severely compromised by 
SOX9 loss. In wild-type HFs, SOX9 marks the HF-SCs 
in the bulge and the upper ORS but is down-regulated as 
HF lineages progress. In contrast, SOX4 and SOX6 are 
expressed in the HG, ORS, and matrix (Rendl et al. 2005; 
Greco et al. 2009; Lien et al. 201 1), and SOX18 and SOX21 
are expressed later in hair differentiation (Pennisi et al. 
2000; Kiso et al. 2009). Whether these changes in gene 
expression compensate and contribute to the relative 
normalcy of Sox9-null HF-SC progeny is beyond the 
scope of the present study. However, in some other 
tissues, redundancy between SOX9 and some of its family 
members has been described (Lefebvie et al. 1998; Stolt 
et al. 2003). 

Adult SCs reside in a discrete microenvironment, 
which plays a dominant role in regulating SC mainte- 
nance, self-renewal, activation, and differentiation. In the 
HF-SC niche, signals from the inner bulge layer, DP, 
dermal fibroblasts, and adipocyte precursor cells have 
all been shown to impact HF-SC quiescence and hair 
cycle activation (Botchkarev et al. 1999; Plikus et al. 
2008; Rendl et al. 2008; Festa et al. 2011; Hsu et al. 2011). 
While these collective studies provide insights into how 
sternness is influenced by extrinsic signaling, few studies 
have focused on whether SCs themselves might express 
factors that impact their own and their neighbors' behav- 
ior and, if so, how. 

After identifying genes whose expression had changed 
upon conditional ablation of Sox9 in adult HF-SCs, we 
realized that a remarkable number of these genes encode 
either secreted proteins or proteins with an extracellular 
domain. Moreover, in HFs containing a mixture of wild- 
type and Sox9 mutant HF-SCs, IFE differentiation was 
induced by not only the mutant but also the wild-type 
HF-SCs. These findings expose a hitherto unrecognized 
communication among HF-SCs within their native niche 
to maintain their own stemness. Moreover, our findings 
also illuminate the importance of SOX9 to niche micro- 
environment. 

It was notable that many of SOX9's bona fide target 
genes encode a group of proteins known to enhance 
Activin/TGFp signaling and that Activin B injections 
rescued the IFE differentiation defect in the Sox9 mutant 
niche. These findings underscore Inhbb, which encodes 
Activin B, as a physiologically significant SOX9 target. 
Moreover, given the striking differences between condi- 
tional loss of TGFp and Activin receptor signaling with 
regard to their effects on the HF-SC niche, our findings 
further suggest that despite their shared ability to acti- 
vate pSMAD2-mediated transcription, TGFps and Ac- 
tivins have nonoverlapping functions that cannot be com- 
pensated for by the expression of each other. In this 
regard, it is intriguing that TGFp2 is induced by HG at 



anagen onset, enhances pSMAD2 in the HG, and pro- 
motes HG proliferation by lowering BMP signaling 
(Oshimori and Fuchs 2012). In contrast, consistent with 
the conventional role of TGFp family members in 
restricting keratinocyte proliferation, Sox9 deficiency 
(and reduced Activin/pSMAD2) was manifested in the 
activated bulge HF-SCs, which did not return to quies- 
cence in mid-anagen. 

Broader mechanistic issues for regalatiag SC behavior 

The complexities that we unearthed in examining one 
SOX9 target, Inhbb, in this study make it unlikely that 
a few linear pathways will explain SOX9's governance of 
stemness. Further support for this notion stems from the 
fact that even though Lhx2 is one of the genes most 
markedly repressed upon loss of SOX9, its loss-of-func- 
tion phenotype is quite distinct from that of Sox9. On the 
one hand, LHX2 and SOX9 orchestrate quiescence and 
stemness, and their loss leads to proliferation and loss of 
HF-SC markers and behavior. On the other hand, each has 
their own impact on lineage suppression, with LHX2 
playing a more important role in sebocyte differentiation 
(Folgueras et al. 2013), and SOX9 playing a greater role in 
suppressing IFE differentiation. Together, these factors 
maintain stemness and direct HF-SCs along the correct 
lineage. That said, our studies show that SOX9 binds and 
regulates Lhx2 in HF-SCs, while LHX2 ChlP-seq did not 
reveal binding to Sox9 (Folgueras et al. 2013). This could 
account for the more severe, combined Activin B/LHX2- 
deficient-like phenotype of our Sox9 cKO mice. 

In closing, we are left with the view that adult HF-SC 
behavior is determined by inputs from multiple niche cell 
types, including SCs themselves. This guarantees that 
the decision of when to make tissue and what lineage to 
follow happens in synchrony and only when a variety of 
cellular inputs from the tissue are met. Our findings 
might also be relevant to human epidermal cells, which, 
unlike their mouse counterparts, express SOX9, which, 
when overexpressed, causes reduced expression of differ- 
entiation-related genes (Shi et al. 2013). 

Materials and methods 

Mice and labeling experiments 

Sox9f"" mice (Akiyama et al. 2002) were crossed with transgenic 
KitlS-CrePGR mice (Morris et al. 2004) and R26YFP mice 
(Srinivas et al. 2001). Cre-PGR was activated in 2-mo-old mice 
by RU486 treatment: 7-d intraperitoneal injection of 1 % RU486 
(VWR) in sesame oil (0.75mg/10g body weight per day; Sigma) 
together with 2-wk topical application of 4% RU486 in ethanol. 
Knock-in Sox9-creERT2 mice (Soeda et al. 2010) were crossed 
with R26YPP mice. CreER was activated by intraperitoneal 
injection of 2-mo-old mice with 2% tamoxifen (Sigma) in corn 
oil (lmg/10 g body weight per day for 2 d; Sigma). For EdU pulse 
experiments, mice received intraperitoneal injections of EdU (25 
jjLg/g; Invitrogen) for 24 h (anagen HFs) or 10 d (telogen HFs). 
Depilation of mid-dorsal HFs was achieved on anesthetized mice 
with molten wax, which was peeled off after hardening. 
Recombinant human Activin B (500 ng; R&D Systems) or BSA 
control was repeatedly intradermally injected for 3 d with Fluo- 
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Spheres (Invitrogen) into the back skin post-depilation. Postnatal 
day 50-60 (P50-P60) CDl mice from Charles River Laboratories 
were used for SOX9 ChlP-seq and ChlP-qPCR experiments. Mice 
were housed and bred in the Association for Assessment and 
Accreditation of Laboratory Animal Care (AAALAC)-accredited 
Comparative Bioscience Center (CBC) at The Rockefeller Uni- 
versity in accordance with National Institutes of Health guide- 
lines, and procedures were performed with Institutional Animal 
Care and Use Committee (lACUC)-approved protocols. 

ChlP-seq assay and analysis 

All materials, methods, and sequencing for ChlP-seq have been 
previously described (Lien et al. 201 1 ). Independent immimopre- 
cipitations were performed on FACS-sorted populations. Iso- 
lation of HF-SCs and IFE cells by FACS has been described 
elsewhere (Nowak and Fuchs 2009; Folgueras et al. 2013). Two 
replicates of SOX9 ChlP-seq were performed in HF-SCs, and one 
replicate was performed in IFE cells; 5 X 10'' cells were used 
for each ChlP-seq run. Rabbit polyclonal anti-mouse SOX9 
(Millipore, AB5535) was used for ChlP-seq. ChlP-seq reads were 
ahgned to the mouse genome (mm9, build 37), and peaks were 
called with MACS software (Zhang et al. 2008). SOX9-bound 
targets were defined as mouse RefSeq genes with a peak at their 
promoters (-2 kb to +2 kb of transcription start sites) or 
enhancers (-50 kb of transcription start sites, gene body, -h5 kb 
downstream from the gene). For the identification of SOX9 
targets, only peaks that were specific to HF-SCs were called. 
The MEME software suite (Bailey et al. 2009) was applied to the 
sequences under the SOX9 peaks (200 bp around summit) for 
finding enriched motifs. ChlP-seq raw data will be submitted to 
Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo). 

Statistics 

To determine the significance between two groups, comparisons 
were made using impaired two-tailed Student's t-test in PrismS 
(GraphPad software). For all statistical tests, the 0.05 level of 
confidence was accepted for statistical significance. 
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